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ABSTRACT 

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in 
the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S'*MC). We detect the 
mid-IR emission from poly cyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly 
increasing the range of environments where PAHs have been spectroscopically detected in the SMC. 
We investigate the variations of the mid-IR bands in each region and compare our results to studies 
of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH 
emission in the SMC is characterized by low ratios of the 6—9 /im features relative to the 11.3 fim 
feature and weak 8.6 and 17.0 /^m features. Interpreting these band ratios in the light of laboratory 
and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those 
in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution 
shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead 
reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of 
our observations for our understanding of the PAH life-cycle in low-metallicity galaxies — namely that 
the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and 
therefore being more susceptible to destruction under typical interstellar medium conditions. 
Subject headings: dust, extinction — infrared:, ISM — Magellanic Clouds 



1. INTRODUCTION 

The mid-infrared emission from normal star-forming 
galaxies is characterized by a series of bright emis- 
sion bands at approximately 3.3, 6.2, 7.7, 8.6 and 11.3 
fim, with weaker features at surrounding wavelengths. 
Laboratory work has shown a close correspondence be- 
tween the interstellar emission features and the vibra- 
tional modes of the carbon-carbon (C-C) and carbon- 
hydrogen (C-H) bonds in polyc yclic aromatic hydrocar- 
bons (PAH s; |Leger fc Puget|[l984[ |Allamandola et ah) 
19851 11989 ). Although there is not a definitive labora- 



tory identification of the mid-IR band carrier or popu- 
lation of carriers, it is considered very likely that these 
features arise from the vibrational de-excitation of PAHs 
with a range of sizes after ab sorbing UV o r optical pho- 
tons (for a recent review, see Tielens|[2008 ). 

The widespread observation of the emission bands in 
the Milky Way and elsewhere suggests that PAHs are 
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an abundant and energetically important component of 
interstellar dust (emitti ng ^ 10% of the total infrared 
emissio n from a galaxy; |Helou et al.||2000] [Smith et al.| 



2007b I . In addition, PAHs play a number of important 
roles in the interstellar medium ( ISM). They are a major 
source of photoelec t ric heating ([Bakes fc T ielensI |1994{ 
I Wolfire et aT][T995l [Hollenbach fc I'ielen s 1999; Hclou| 
let al.l2001 [) and participate in chemical reaction networks 



(Bakes & Tielens 1998 



Wakelam & 



^ [Wolfire et al.||2008[ 

|Herbst(2005 ) in a variety of ISM phases. Because of these 
crucial roles, we would like to understand how and why 
the characteristics of PAHs depend on galaxy properties 
such as metallicity and star-formation history and what 
specific factors control the abundance, size distribution 
and physical state (ionization, hydrogenation, compact- 
ness, symmetry, etc.) of PAHs in the ISM. 

The physical state of a population of PAHs is reflected 
in the relative intensities, profile shape and central wave- 
length of the mid-IR emission bands. Of the major PAH 
bands, laboratory spectroscopy has suggested that those 
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at 3.3, 8.6 and 11—14 /im arise from C-H stretching and 
bending modes while the 6.2 and 7.7 complexes are from 
C-C stretching. Longer wavelength features like the 17.0 
/xm complex are most likely rel ated to more comple x ring 
bending modes of the molecule ( Moutou et al. 2000 1 . The 
relative strengths of these bands tor a given population 
of PAHs will depend on their size distribution and phys- 
ical state allowing us diagnostics of the PAH population 
from mid-IR spectroscopy. 

For instance, laboratory and theoretical investigations 
have consistently shown that the C-H out-of-plane bend- 
ing modes between 11—15 /zm in ionized PAHs are weak, 
while they are strong in neutral PAHs. The opposite is 
true of the C-C stret ch and C-H in-plane bending modes 



et al. 



between 6 — 9 /zm (Szczepanski & Vala 1993 Hudgins 



1994|, so the ratios ot these two sets of features 
to some degree. Recent 
argued that the varia- 



should trace PAH ionization 



work by Galliano et al 
tions of the 6.2, 7.7, 



( |2008a 

^6 and 11.3 /zm features within 



and between galaxies are consistent with being controlled 
by ionization. Specifically, they find that the ratios of 
6.2, 7.7 and 8.6 to the 11.3 /xm bands are correlated 
over an order-of-magnitude, while ratios such as 7.7/6.2 
and 8.6/6.2 stay roughly constant as a function of the 
7.7/11.3 ratio. 

The ratios of short-to-long wavelength PAH bands, 
however, are also expected to d epend on the PAH size 
distribution (Draine & Li 2001), such that small PAHs 
emit more strongly in shorter wavelength features. In 
addition to this general trend, laboratory and theoretical 
studies have suggested t hat some specific fea tures orig- 



inate from large PAHs. Peeters et al. (2004) find that 



large PAHs are required to reproduce the general char 
acteristics of the 17.0 /im complex, although the charge 
sta te of the carrier is not wel l deter mined. Recent work 
by Bauschlicher et al. (2008 2009) has suggested that 



the strength ot the 8.6 /xm band also is an indicator 
of the contribution of large PAHs to the mid-IR spec- 
trum. They find that in their theoretical investigations 
of the spectrum of large compact and irregular PAHs, 
only large, charged, compact PAHs have a feature at 8.6 
/im that can reproduce the interstell ar band. 



Fi nally, as reviewed i n det ail by Hony et al. (20011 
and [Bauschlicher et al. (2009), the relative streiigths of 
the 11—15 /im C-H out-ot-plane bending modes provide 
information about the structure of PAHs. The type of 
structure adjacent to the vibrating C-H bond determines 
the wavelength of the mode. If a given aromatic carbon 
ring on the edge of a PAH has only one attached hydrogen 
(out of four possible for a PAH that consists of more than 
one aromatic ring) , while the rest of the bonds arc con- 
nected to other aromatic rings, it is considered a "solo" 
hydrogen. If two of the possible sites contain C-H bonds, 
this is considered a "duo" , and the "trio" and "quartet" 
are similarly defined. In structural terms, straight edges 
produce "solo" hydrogens while corners produce "duo" 
hydrogens. Irregularities like bays and linear extensions 
produced "trio" and "quartet" hydrogens. These differ- 
ent modes have diff erent wavelengths (depen dent on the 
size of the molecule; Bauschlicher et al.|2009 ). Thus, ob- 
servations of the strengths of the 11—15 /xm features can 
provide information about the PAH structure. 

Using the knowledge of PAH bands gained from these 
laboratory and theoretical studies, it is possible to di- 



agnose the physical state of PAHs in extragalactic envi- 
ronments observed with ISO and Spitzer. In particular, 
evidence for changes in PAH properties with radiation 
field intensity and metallicity have been seen in a variety 
of studies: 



• Galliano et al. (2008a I studied the variation of 
the ratios of the 6.2, 7.7, 8.6 and 11.3 bands in 
a diverse sample of galaxies and individual star- 
forming regions and found variations in the 6—9 
vs 11.3 /im features correlated with the local UV 
field strength — a trend they interpreted as arising 
mainly from changes in the PAH ionization state. 



Smith et al. (2007b) have found that PAH emis- 
sion from A(3M -hosting galaxy nuclei shows very 



low 7.7/11.3 ratios (see a lso O'Dowd et aT] 2009 



Rieke 2010^ Wu et al.||2010p. 



Diamond-Stanic &: ^ ^ ^ 

Some studies have also found that the spectra of 
these sources also show average or even high ra- 
tios of the 17.0 feature to the 7.7 and 11 .3 bands 
(iSmith et al.||2007b| lO'Dowd et al. 2009). These 



trends have been interpreted as the preferential de- 
struction of small PAHs by the hard radiation field 
produced by the AGN, leading to decreased emis- 
sion in the short wavelength bands relative to the 
long wavelength bands. 



• Hunt et al. (2010) studied the PAH features in 
a sample of blue compact dwarfs (BCDs) with a 
range of metallicities and found evidence for in- 
creased strengths of the 8.6 and 11.3 features rela- 
tive to the other PAH bands in these intensely star- 
forming galaxies. They interpreted these trends 
as evidence for destruction of small PAHs in the 
harder, more intense radiation fields of these galax- 
ies cit i ng re cent theoretical work by |Bauschlicher| 



et al. (2008) which suggests that the 8.6 /xm fea- 



ture arises predominantly from large PAHs. 



• Smith et al. (2007b) measured the integrated PAH 
band strengths in the central regions of galaxies in 
the SINGS sample and found indications that the 
17.0 /im complex is weak compared to the other 
PAH bands at low-metallicity, although their sam- 
ple contains few low-metallicity galaxies. 

• Resolved studies of PAH emission in the Magellanic 
Clouds have also sh own variations in the band ra- 
tios. For instance. Reach et al. ( |2000 ) detected 
PAH emission from the quiescent nioiecular clouds 
SMC Bl #1 with a very low 7.7/11.3 ratio. |Li 
fc Draine| |2002 ) found that they could not repro- 
j 4.u,„ uit;i^ 1 ratio using their Milky Way PAH 



duce the S 
model 



In addition to these changes in the PAH physical state, 
observations with ISO and Spitzer have demonstrated 
that the fraction of PAHs relative to the total amount 



of dust decreases in low-metallicity environments (|Mad 
denl|2000l [Engelbracht et al.||2005l [Madden et al. 



Wu et al.|r2006; Draine et al. 



120071 



Hunt et al.[[2010[ iLeboutciller et a. 



Gordon et al 
201ip . A nu"^ 



2008 



ber o ' 

explanatioiis tor tlEe PAH deficit in low-metallicity galax- 
ies have been proposed including a time lag between the 
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formation of PAHs in AGB star atmosph eres relative to 
SNe-produced dust ( Galliano et al. 2008b|), destructi on of 



PAHs by supernova sliocks ( |U'HaIloran et a"L||2006[ ) and 
destru ction of PAHs by ha rder and/ or more intense UV 



fields (Madden et al. 2006; Gordon et al. 20081. Recent 



work by Sandstrom ct al. (2010i) (Paper 1) on the frac- 
tion of PAHs in the Small Magellanic Cloud has shown 
a very low PAH fraction in the diffuse ISM of this low- 
metallicity galaxy and a PAH fraction in dense gas within 
a factor 2 of the Milky Way. They argue that PAHs 
are destroyed in the diffuse ISM of the SMC and may be 
forming in dense regions. Thus, in addition to variations 
in the destruction of PAHs, it is possible that as a func- 
tion of metallicity the balance between AGB-produced 
PAHs and PAHs that form in the ISM may change. 

These proposed variations in the life-cycle of PAHs in 
low-metallicity galaxies will be reflected in changes to the 
physical state of the PAHs themselves — their size dis- 
tribution, hydrogenation, ionization state, chemical sub- 
stitutions and/or functional groups and structural char- 
acteristics. Enhanced destruction of PAHs by UV fields 
or shocks should change the PAH size distributio n, since 
smaller molecules are easier to phot odissociate ( Allain 



et al. 1996a|b Le Page et al. 2001). If the dommant 



formation mechanism of PAHs changes as a function of 
metallicity it may be refiected in the size distribution 
and irregularity of the molecules. Thus, by observing 
the physical state of PAHs in a low-metallicity environ- 
ment we hope to understand what drives the observed 
differences in the PAH life-cycle. 

To that end, we present the results of a mid-IR spec- 
troscopic survey of PAH emission in the Smal l Magel- 
lanic Clou d (SMC). At a distance of 61 kpc ([Hilditch 
etal.|2005D and with a metallicity of 12+ log(0/Hj ~ 8.0 



( Kurt fc DufourTr998[ ) , the SMC presents the ideal target 
to understand the properties of PAHs in low-metallicity 
systems at high spatial resolution. The SMC represents 
an intermediate case between Galactic sources where 



varia t ions can be studied within objects ( Hony et al. 
2001] JPeeters et al.|[2002l [Rapacioh erar]|2005| igmle 
et al. 2007[ ) and other galaxies where only galaxy-wide 
results are obtained. The Spitzer Spectroscopic Survey 
of the SMC (S'^MC) consists of low-resolution spectral 
mapping observations using the Infrared Spectrograph 
(IRS) on Spitzer over six of the major star-forming re- 
gions of the SMC. In Paper I, we presented the results of 
SED and spectral fitting using the S'^MC data to deter- 
mine the PAH fraction across the galaxy. In this paper 
we turn our attention to using the band ratios to learn 
about the physical state of PAHs. 

In Section [2] we describe the S^MC observations and 
data reduction. Section |2.5| covers the spectral fitting. 
In Section [3] we describe the variations in the band ratios 
in the six star-forming regions in our sample and finally 
in Section |4] we discuss what the variations tell us ab out 

the 
PAH 



the physical state of SMC PAHs and in Section |5.1 
implications of our results for understanding the 
deficiency in low-metallicity galaxies. 



2. SPITZER SPECTROSCOPIC SURVEY OF THE 
SMALL MAGELLANIC CLOUD (S-^MC) 

We mapped six star-forming regions in the SMC us- 
ing the low-resolution orders of the InfraRed Spectro- 
graph (IRS) on Spitzer as part of the S'^MC project (GO 



30491 — now combined with a larger mosaic of the galaxy 
as part of SAGE-SMC; |Bolatto et al.|20071 Gordon et al. 
2011, submitted). Figure [I] shows the coverage of our 
spectral maps ov erlaid on tne S'^MC 24 /im mosaic (Bo 



latto et al 
images of t 



20071 



Figures [2] [3] and H show three-color 
tie SIX regions, overlaid with the coverage of 
all of the orders. These regions are all displayed with 
the same color scale. In the following analysis, we are 
limited to the regions with coverage in both of the SL 
orders because of our focus on the PAH emission bands. 
We list the details of the observations in Table [T] 

The regions we have mapped were selected to cover 
a range of properties and therefore represent a diverse 
sample of SMC sources. The N 66 field covers the well 
studie d giant H II region NGC 346, which hosts 33 O 
sta rs (Klassey et al. l 1989 ). ISO observations of N 66 
by Contursi et ai. ~( |2boo[ ) show radiation field intensity 
upwards of 10^ times the solar neighborhood radiation 
field in some of the brightest mid-IR peaks. We have 
also mapped the quiesc ent molecular regio n SMC Bl #1, 



found in the SW Bar ( |Rubio et al.||1993i ). This was the 
first location where PAh emission was dete cted spectro- 
scopically in the SMC bylReach et al.|(|2000|). Also in the 



SW Bar region, we map three other fields: N22, a com- 
pact H II region surrounded by molecular gas; SW Bar 
1, a region of active, but more evolved, star-formation 
(including the N 27 nebula), with bright 8 /im and Ha 
emission; and SW Bar 3, which covers regions of molec- 
ular gas and recent star-formation. In the Wing of the 
SMC, we map the N 83/N 84 complex, which covers a re- 
gion oniigh_UVfieldintensity with abundant molecular 



gas (Leroy et al. 2009). Finally, in the N 76 region our 



g i 

( Stanimirovic et al.|2005] [Sandstrom et al |2009[ ), a Wolf- 
Rayet nebula and a few small regions of molecular gas 
on the outskirts of the nebula. 

2.1. Mapping Strategy and Spectral Cube Construction 

The maps were constructed by stepping the IRS slit 
in a grid of positions parallel and perpendicular to the 
slit. The maps are fully sampled by stepping by half of 
a slit width (steps of 1.85" and 5.08" for SL and LL, 
respectively) perpendicular to the slit. The LL maps 
are composed of a grid of 98 perpendicular and 7 par- 
allel steps, with half slit length steps parallel to the slit 
and 14 seconds of integration time per position (except 
N 76 which has 75 x 6 steps perpendicular and paral- 
lel). The LL spectral maps cover an area of 493" x 474" 
(376" X 395" for N 76). The SL maps are constructed 
similarly except that we use full slit width steps parallel 
to the slit to increase the mapped area at the expense of 
some pixel redundancy. All of the SL maps, aside from 
the map of SMC Bl, have 120 perpendicular by 5 par- 
allel pointings (220" x 208") with integration times of 
14 seconds per position. We have made a deeper map 
of SMC Bl, the location where PAH emission was first 
detected from the SMC fieach et ar '2000) using 60 by 
4 slit positions (109" x 156") with 60 second integration 
times. The spectral coverage of the low-resolution orders 
extends between 5.2 and 38.5 /im with spectral resolving 
power ranging between ~ 60 — 120. 

Each mapping observation was preceded or followed 
closely in time by an observation of an "off" position 
at RA lh9"M0" Dec -73°31'30" (J2000). This posi- 
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Fig. 1. — The coverage of the S'^MC survey overlayed on the MIPS 24 fim map from S^MC. The color scale is logarithmic, with the 
stretch illustrated in the colorbar. The red boxes show the coverage of the LLl order maps (the LL2 maps are shifted by ^ 3') and the 
green boxes show the coverage of the SLl order maps (the SL2 maps arc shifted by ~ 1'). We also identify the various regions of the galaxy 
by the names we will refer to in the remainder of this paper. 



TABLE 1 

S^MC Observations 



Region R. A. (J2000) Dec. (J2000) Observation Date AOR Number 



LL Observations 



SW Bar 3 


0^45^23= 


-73°22'00" 


11 Jun 2007 


18264832 


N 22 


Qh47m42s 


-73°15'06" 


09 Sep 2006 


18261504 


SW Bar 1 


0''48™08= 


-73°08'24" 


11 Jun 2007 


18264576 


N 66 


0''59™09= 


-72°09'14" 


15 Nov 2006 


18262272 


N 76 


1^03^41= 


-72°01'30" 


12 Dec 2006 


18263040 


N 83/84 


lhl4™05= 


-73° 17' 00" 


10 Sep 2006 


18263808 


SL Observations 


SW Bar 3 


0^45^26= 


-73°22'30" 


10 Jun 2007 


18265344 


SMC Bl #1 


0''45™21= 


-73°20'10" 


10 Jun 2007 


18261248 


SW Bar 1 


0''48™13= 


-73°05'20" 


10 Jun 2007 


18265088 


N 22 North 


0''48™11= 


-73°13'56" 


20 Nov 2006 


18262016 


N 22 South 


0iM8™04= 


-73°17'50" 


20 Nov 2006 


18261760 


N 66 North 


0'>59™22= 


-72°09'07" 


21 Nov 2006 


18262784 


N 66 South 


0'>59™05= 


-72°12'50" 


20 Nov 2006 


18262528 


N 76 1 


l''03™ll'= 


-72°04'01" 


09 Dec 2006 


18263296 


N 76 2 


l''04™00= 


-72°02'50" 


09 Dec 2006 


18263552 


N 83 


1^14^00= 


-73°16'40" 


10 Dec 2006 


18264064 


N 84 


1^14^50= 


-73°19'20" 


10 Dec 2006 


18264320 



SMC PAH Spectroscopy 



5 




05'0' 



lO'O" 



-72°15'0' 




oroo^ao' 



00"59"30' 



0= 




Fig. 2. — Coverage of the S^MC observations in the N 22 and N 66 regions. The three-color images are constructed from S^MC 
photometry in the IRAC 3.6 (blue), IRAC 8.0 (green) and MIPS 24 (red) bands. The color table is the same in Figures[2] [sjand [i] The 
coverage of the LLl and SLl orders is shown with a solid white line while the LL2 and SL2 coverage is shown with a gray dashed line. 
For the following analysis we are limited to regions with coverage in both SL bands. The white boxes illustrate locations where we have 
extracted spectra from the cube for this analysis (see the text for details) and the right panel shows all of the extracted regions along with 
a number identifying them. There is a very bright point source in the N 22 region which is saturated in the 24 fim MIPS observation, we 
show its location with a black star. Because of the artifacts associated with bright point sources we are forced to omit some regions from 
our analysis. 
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-71°55'0' 



OO'O' 



05'0' 



-72°10'0' 





oros^o^ 



oa^o' 



lO'O' 



15'0" 



20'0" 



-73°a5'0' 





0ri6"0n5"30 



Fig. 3. — Coverage of the S*MC observations in N 76 and N 83. The images and color table are the same as in Figure [S] In N 76, we 
avoid the region dominated by the supernova remnant IE 0102.2— 7219 (visible as a ring-like feature in the 24^m band; Sandstrom et al.| 
|2009[ presented an analysis of the IRS observations for this region). In N 83, most of the area of the northern SL cubes are atlected by 
artifacts due to the saturation of the peak-up array due to a very bright point source in the region. We avoid those areas in selecting the 
extraction regions. 



tion was chosen to be free from SMC emission at the 
level of our MIPS and IRAC mosaics from S^MC. The 
"off" position was observed with the same ramp times 
as its corresponding mapping observations using "star- 
ing mode". There were four nod positions in each "off" 
observation, with 9 repetitions per nod. In the course 
of processing the maps, the outlier-clipped mean of the 
36 two-dimensional spectral images of the "off" position 
was subtracted from each of the mapping BCDs. This 



removes the zeroth-order foreground emission from the 
zodiacal light and Milky Way cirrus as well as mitigat- 
ing the effects of time-variable rogue pixels in the low- 
resolution arrays. 

We have used the BCDs from pipeline version S15.3 
and S16.1 to create our final cubes. There are no sig- 
nificant differences between these two pipeline versions 
for our purposes. The cubes were assembled using the 
Cubism software package (vl.7) distributed by the SSC 




15'0' 



20'0" 



-73°25'0' 




SMC Bl #1 





Bar 3 



00"47 



Fig. 4. — Coverage of the S'*MC observations in the SW Bar 1 and 3 regions. The images and color table are the same as in Figure [2] 
The coverage of the SMC Bl spectral map is shown in the SW Bar 3 region. 



(Smith et al. 2007a). Cubism takes the pipehne pro- 
cessed BCDs as mput and uses polygon- cUpping based re- 
projection to transform the two dimensional spectral im- 
ages into a three dimensional spectral cube. Cubism also 
propagates the uncertainties associated with the data 
collection events, known as ramps. We propagate these 
uncertainties throughout the following analysis, but note 
that they represent only the statistical uncertainties of 
each pixel and do not take into account the more im- 
portant systematic errors. In assembling the cubes in 
Cubism we apply a "slit-loss correction function" which 
is essentially an extended source calibration factor for 



the IRS dSmith et al.|[2007a |. 

We reniove bad pixels in tne cubes using the tools pro- 
vided in Cubism. These bad pixels are identified by the 
stripes they produce in the spectral maps, resulting from 
the movement of the bad pixel through the grid of po- 
sitions. The number of bad pixels increases towards the 
longest wavelengths in the cube, effectively making the 
spectrum past ~ 35 fim have much lower signal-to-noise. 
Regarding bad pixel flagging, we note that in cases where 
the core of the PSF of a point source falls on a bad pixel 
in the detector, there can be difficulties in reconstructing 
the emission profile in Cubism. This issue is exacerbated 
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by two problems: the spatial under-sampling of the IRS, 
and the "full-width" steps used to construct the map, 
which reduces the pixel redundancy. Because we are in- 
terested in the diffuse emission, we have performed ex- 
tensive bad pixel subtraction in order to get the cleanest 
map of the region without regard to the effects on point 
sources. For that reason, in the vicinity of bright point 
sources there can be artifacts that interfere with the de- 
termination of the PAH emission band strengths in the 
spectra of point sources . In the following, we avoid these 
regions in the cubes, particularly in the core of N 66 and 
in the vicinity of the bright point sources in N 22 and N 
83. In N 83, the brightest point source presents another 
difficulty in that it saturates the peak-up array, leading 
to uncorrected droop effects over a large fraction of the 
map. In our analysis, we have avoided regions affected 
by this artifact as well. 



2.2. Foreground Subtraction 

The "off" position observations that were performed 
close in time to each mapping BCD allow us to subtract 
the majority of the zodiacal and Milky Way foregrounds 
from our maps. There are, however, gradients in both 
foregrounds between the "off" location and the map lo- 
cation. To remove these additional gradients we use the 
zodiacal light spectrum predicted by the Spitzer Observa- 
tion Planning Tool (SPOT) at the times of each observa- 
tion to calculate the difference in the zodiacal foreground 
between the "off" location and the map location versus 
wavelength. SP OT uses the zodiacal light model of Kel- 
sall et al. (1998) from DIRBE observations. The modifi- 
cations to this model for predicting the zodiacal light at 
Spitzer wavelengths are described in the documentation 
from the Spitzer Science Center(SSC)^. To determine the 
gradient in the Milky Way foreground we assume a cor- 
relation between the cirrus dust e mission and the column 
density of n eutral hydrogen (as in Boulanger et al.|1996 ) . 
We use the Draine & Li (2007) emissivity of Milky Wa,y 
dust and a map of H I from a combined ATCA/Parkes 
survey (E. MuUer, private communication) to measure 
the Milky Way foreground gradient as in Paper I. 



2.3. "Dark Settle" Related Artifacts in the S^MC 
Spectral Mapping 

In performing the analysis of the S'^MC data we discov- 
ered that both the SL and LL observations are affected 
by an artifact related to a time- and position variable 
residual background level on the detector. This artifact 
is a general issue in IRS spectral mapping, but is particu- 
larly important for faint regions like those we study here. 
This issue is very similar to what has been seen in LH 
observations and referred to as the "dark settle" issue^. 
The artifact causes a tilt to the various orders, leading 
to a mis- match between the spectral orders LL1-LL2 and 
SL1-SL2 as well as introducing spurious curvature to the 
spectrum. Because of the faintness of the PAH emission 
in many SMC regions, this artifact can have a significant 
effect on the flux in the 7.7 /im complex, located at the 
overlap of the SLl and SL2 orders. 
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features? 



Our basic approach to dealing with the artifact is to 
correct the data at the BCD level by determining the 
shape of the residual background from the inter-order 
region. The details of this procedure are described in 
the Appendix. For the SL orders in particular, the inter- 
order space is very small, especially since the gap between 
the SL3 and SL2 orders is strongly affected by the scat- 
tered light from the peak-up arrays and its correction in 
the IRS pipeline. At a basic level, we do not have infor- 
mation about the background level within the spectral 
orders, which is precisely where we would like to remove 
it. Therefore, our correction of this artifact is approx- 
imate at best and we proceed in the remainder of the 
paper to deal with two versions of the spectral cubes: 
one with no correction of the "dark settle" issue and one 
with a correction we determine from the inter-order light. 
As described in the Appendix, the technique we use may 
overcorrect (i.e. by removing too much flux from the 
7.7 /im feature). We argue that the uncorrected and cor- 
rected spectra should bracket the range of possible values 
for the band strengths. 

2.4. Spectral Extraction and Stitching 

In order to get better signal-to-noise for determining 
the band ratios we extract the spectra in 25" x 25" boxes. 
This extraction aperture is also large enough to avoid 
problems with the variation of the PSF with wavelength 
(the FWHM of the PSF at 35 pcm is ~ 9"). We have 
tiled the cubes in the region of full overlap between all 
orders with these extraction boxes. The right panels of 
Figures [2j [3] and H] show the extraction boxes for each 
region, excluding those affected by point source artifacts. 
We refer to the individual spectra by the number shown 
on these diagrams throughout the paper. 

After the correction for foreground gradients and "dark 
settle"-type artifacts, our extracted spectra show offsets 
in the overlap of the SLl and LL2 orders that vary from 
region to region, but stay relatively constant within a 
given region. These offsets are additive and on the order 
of a few tenths of MJy sr~^. In correcting these offsets 
we first extract photometry in matching apertures from 
the S'^MC 24 /im mosaic (wh ich has been f oreground sub- 
tracted as described in jSandstrom et al. 2010) and find 
the average offset at 24 /tm for the MIPS and synthetic 
IRS photometry for each region. We then use this ad- 
ditive region-based offset to tie the IRS spectrum to the 
MIPS photometry — a step which is necessary because we 
have no prior information on which order is "correct" and 
the 8 pLTci photometry relies on an extended source cor- 
rection for IRAC and is thus less reliable than the 24 /im 
photometry. Next we determine the offset between the 
SL and LL orders in their overlap region (~ 14.3 — 14.7 
/tm) and find another additive region-based correction 
that we then add to the SL spectra. 

After tying the spectra to the MIPS 24 /tm photometry 
and correcting the LL/SL offset, we stitch the spectra to- 
gether by interpolating the long wavelength end of each 
order onto the grid of the overlapping shorter wavelength 
order and averaging the resulting values together. The 
procedures we employ in correcting offsets between the 
orders and stitching the spectra together do not signifi- 
cantly affect our results. 

2.5. Fitting the PAH Emission Bands 
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We use th e PAHFIT spectral fitting routine (Smith 
et al.||2007b[ ) to measure the strengths of the PAH emis- 
sion bands and a variety of emission Hnes from our spec- 
tra. Figure [5] shows the resuhs of fits to some high signal- 
to-noise spectra from our dataset. PAHFIT uses combi- 
nations of Drude profiles to fit the emission bands and 
Gaussians to fit the e mission lines. We do not include 
extinction in our fits — Lee et al. (2009) found very low 
levels of mid-IR extinction in the SMC as expected given 
the galaxy's low dust-to-gas ratio. The effect of extinc- 
tion on the band ratios would tend to primarily suppress 
emission in the 8.6 and 11.3 /xm bands because of the 
silicate extinction feature at 9.7 fiui. In the cases where 
the PAH emission features form a complex (for instance, 
the 7.7 fim feature), we report only the total and not 
the properties of the individual components. We have 
used the default line list provided in PAHFIT with the 
exception of adding in the Pfund-a recombination line of 
hydrogen at 7.46 /im, which we see in the spectra of N 
66 and some of the spectra in N 22. We perform the fits 
on the original spectra and on the "dark settle" artifact 
corrected spectra. 

Table [2] lists the measured integrated intensities for the 
major PAH features as well as the lines of [Ne HI] at 15.6 
/im and [Ne II] at 12.8 /im for each spectrum extracted 
from the cubes both with and without the "dark settle" 
correction. 

3. RESULTS 

3.1. Band Ratios in the SMC 

In the following section we present our results in a se- 
ries of plots showing the band ratios determined in each 
of the extraction regions shown on Figures [2] |3] and |4j 
For comparison we also plot band ratios from two liter- 
ature samples: 1) the central regions of galaxies in the 
SINGS samp le without evidence for AGN from |Smith| 
et al. ( 2007b ) and 2) band ratios measured in a sa mple of 



blue compact dwar fs (BCDs) and starbursts from Engel 



bracht et al. (2008). These samples are chosen to repre- 
sent a range of metallicities and star-formation rates and 
to have been measured with the same PAHFIT technique 
to ensure consistent results. The spectroscopic sample of 
SINGS galaxies without AGN represent galaxies that are 
generally more me tal-rich than the S MC. The starburst 
sample from Engclbracht et al. ( 2008 1 , on the other hand, 
span a range of metallicities, from Milky Way metallic- 
ity down to lower than that of the SMC. We note that 
our measurements are on small scales inside individual 



star-forming region s, while the Smith et al. (2007b) and 
Engelbracht et al. (2008) measurements encompass large 
regions of a galaxy, 
caution 



Thus, we use the comparison with 
However, the results of Paper I suggest that 
most of the PAHs in the SMC reside in the vicinity of 
star-forming regions where there is molecular gas, and 
the majority of the PAH emission arises in dense PDRs, 
so in terms of the total PAH emission from the SMC, 
star-forming regions like those we have mapped are likely 
to dominate because of their high luminosities and the 
low PAH fraction elsewhere. 

On the band ratio plots (Figures[6]through 11 ) we show 
a separate panel for each SMC region. To be conserva- 
tive given the uncertainties related to the "dark settle" 
artifact, we use 5a upper limits on the band strength. If 



a given region is a 5(t limit in either the uncorrected 
or artifact-corrected spectrum, we consider it a non- 
detection and use the least stringent limit on the cor- 
rected or uncorrected band ratios on the plot. Our choice 
of 5<T upper limits may introduce a bias into our results 
towards the brighter PAH emission regions. Because of 
the variability in position and time of the artifact, we 
cannot make a simple cut to judge the quality of a spec- 
trum, so we have inspected the spectra by hand and as- 
signed as "best" quality (i.e. not strongly affected by the 
artifact) or "normal" quality. The "best" quality points 
are shown in order to highlight our highest confidence 
band ratios least affected by systematic shifts due to the 
artifact. In these figures we show only the uncorrected 
ratios (aside from the case of limits as described above). 
To illustrate the systematic effects of the artifact on the 
spectra we show an arrow on each panel that connects 
the weighted mean'^ of the uncorrected distribution to 
the weighted mean of the artifact-corrected distribution. 
In determining the weighted mean we have used all spec- 
tra where all of the relevant PAH features are detected 
at > 5a in a given region. For this reason, the centroid of 
the distribution changes slightly from figure to figure due 
to the different combinations of bands shown. In panel 
h) of each ratio plot, we show the SINGS and starburst 
samples and the weighted means of the SMC regions. All 
panels show dashed lines indicating the weighted mean 
of the ratio for the SINGS galaxies without AGN. 

In Table [3] we list the weighted mean band ratio for 
each region along with the uncertainty on the mean. 
Note that in this table the band ratios may differ slightly 
from the values shown on the band ratio figures. This 
is because when calculating the values for the figures 
we have done the weighted mean for only those spec- 
tra where all bands shown in the relevant ratios have 
been detected (e.g. 6.2, 8.6, 7.7 and 11.3 for a plot of 
8.6/6.2 versus 7.7/11.3). For the Table, we include in our 
weighted means all ratios where the two relevant bands 
were detected, so these may differ slightly from those 
shown on the figures. 

3.1.1. The 6.2, 7.7, 8.6 and ll.S^im Bands 

Figuresl6j [Tj [8]and|9]show the flux ratios of the major 
PAH bands at 6.2, 7.7, 8.6 and 11.3 /im. In these Figures 
we plot the other ratios versus the 7.7/11.3 ratio, since 
it generally has the highest signal-to-noise. We can see 
from Figuresl6]through|9]that the 7.7/11.3 is consistently 
low in the SMC compared to the SINGS and starburst 
samples with typical values of ~ 2 — 3 compared to the 
SINGS average of ~ 4 — 5. In fact, almost all of the indi- 
vidual extracted spectra show 7.7/11.3 ratios below the 
SINGS average, reinforcing this conclusion. The ratios of 
the 7.7 and 8.6 features to the 6.2 feature are also distinc- 
tively low in the SMC regions compared to the SINGS 
and starburst samples: Figures [6] shows that the 7.7/6.2 
ratio is low in essentially all SMC regions and Figure [7] 
shows similar results for the 8.6/6.2 ratio. 

In Figurejsjwe show the 6.2/11.3 ratio versus 7.7/11.3. 
Interestingly, the 6.2/11.3 ratio in the SMC is much 
closer to, although still slightly below, the SINGS av- 
erage. Comparing this result with the previous ratios 

Throughout this paper, the weighted mean is defined as the 
mean weighted by the uncertainties as l/cr^ . 
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Fig. 5. — These plots show example PAHFIT results for spectra from our sample. The spectra are labeled with the region and a number 
which corresponds to labels on Figures [2] [s] and |4| Each plot shows the spectrum and statistical uncertainties in gray, the full best fit 
model in black, the PAH contribution to the best nt model in green, the line emission in magenta, the stellar continuum in blue and the 
dust continuum in orange. 



involving the 6.2 band, it appears that the low ratios of 
7.7/11.3 and 7.7/6.2 are mostly due to weakness of the 
7.7 feature and not unusual strengths in the 6.2 or 11.3 
features. We will return to this subject later in the dis- 
cussion. Finally, Figure [9] shows consistently low ratios 
of 8.6/11.3, reinforcing the weakness of the 8.6 feature in 
the SIMC regions. 

To summarize our findings on the 6.2, 7.7, 8.6 and 11.3 
jivn features: the SIVIC regions are consistently lower than 
the SINGS average in the 7.7/11.3, 8.6/11.3, 7.7/6.2 and 
8.6/6.2 ratios. The 6.2/11.3 ratio is only slightly lower 
than the SINGS average, suggesting that the band ratios 
involving the 7.7 and 8.6 /im features are low because 
of the weakness of those bands (conversely, however, it 
could be argued that the 6.2 and 11.3 /im features are 
abnormally strong). We note that the correction of the 
"dark settle" artifact does not strongly influence these 
conclusions and in some cases makes the ratios even more 
extreme compared to the SINGS average. We note that 
if extinction were a major concern, the 11.3 /im feature 



would be suppressed, leading to even more extreme band 
ratios compared with the SINGS sample. 

3.1.2. The ll.O^jLm Band 

In Figures [l0| and [TT] we show the ratios of the 17.0 /im 
band to the 7.7 and 11.3 /xm bands. Figure 10 shows that 



the 17.0/11.3 ratio is consistently low in the SMC regions 
compared to the SINGS mean. Aside from a few outliers, 
all of the galaxies from SINGS and the starburst sam- 
ples which cover the same range of 7.7/11.3 have higher 
ratios of 17.0/11.3. There does not appear to be a clear 
trend in the 17.0/11.3 vs 7.7/11.3 ratio, particularly since 
the starburst sample has galaxies with similar 17.0/11.3 
ratios but 7.7/11.3 ratios spanning more than an order 
of magnitude. Because the charge state of the carrier 
of the 17.0 /im complex is not well-constrained, we also 
show the ratio of 17.0/7.7 in Figure 11 The region av- 



erages of the 17.0/7.7 ratio in the SMC are closer and 
at times above the SINGS average, spanning a range of 
~0.06— 0.1. We note that for these ratios, the artifact 
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Fig. 6. — The 7.7/6.2 ratio versus the 7.7/11.3 ratio. The symbols for each region are indentified in the bottom-right corner of panel g. 
A value is considered a limit in this plot if either the uncorrected or artifact-corrected spectrum does not have a detection of all bands at 
> 5(7. In that case, we use the least stringent of the corrected or uncorrected ratios to be conservative. The "best" and "normal" quality 
points are assigned via visual inspection of the spectra to give some idea of which points are strongly affected by the artifact. The red arrow 
connects the weighted mean of the uncorrected band ratios to the same value for the corrected ratios. These weighted means are shown 
also in panel h) with their colors identified by a representative point in each panel a) through g). The dashed line shows the weighted mean 
values for the SINGS galaxies. Panel h) shows the ratios for the SINGS sample without AGN and the starbursts from [Engelbracht et al.| 
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TABLE 2 
Spectral Fit Results 



Region 


Num 




R.A. 


Doc. 




6.2 






7.7 






8.3 






[Nc II] 




[Nc III] 


N 


22 


1 


11 


.891834 


-73, 


.313254 


4, 


.237±0, 


.527 


5, 


.691±0, 


.778 


0, 


.361±0, 


■ 371 


0^ 


■213it0^ 


,052 


0, 


■ 034it0, 


■ 019 


N 


22 


2 


11 


.905884 


-73, 


.299297 


1, 


,950±0, 


,399 


4, 


.367±1, 


.034 


0, 


■ OOOitO, 


■ 269 


0, 


■040it0^ 


,048 


0, 


■ 003it0, 


■ 022 


N 


22 


3 


11 


.909261 


-73, 


.308280 


2, 


,895±0, 


,455 


3, 


.297±0, 


.939 


0, 


■ OOOitO, 


■379 


0, 


■OSlitO^ 


,075 


0, 


■ 031it0, 


■ 028 


N 


22 


4 


11 


.923294 


-73, 


.294323 


1, 


.998±0, 


.295 


7, 


.586±1, 


.009 


0, 


■ OOOibO, 


■ 331 


0^ 


■222it0^ 


,057 


0, 


■ 046it0, 


■ 020 


N 


22 


5 


11 


.923392 


-73, 


.284345 


6, 


,922±0, 


,409 


14, 


.024±1, 


.014 


1, 


.669±0, 


■ 206 


0, 


■352it0^ 


,057 


0, 


■ 054it0, 


■ 023 


N 


22 


6 


11 


.923489 


-73, 


.274367 


5, 


,214±0, 


,474 


11, 


.883±0, 


.991 


1, 


.311±0, 


■409 


0, 


■301it0^ 


,054 


0, 


■ 069it0, 


■ 021 


N 


22 


7 


11 


.926679 


-73, 


.303306 


2, 


.149±0, 


.448 


5, 


.763±1, 


.168 


0, 


■ OOOibO, 


■ 346 


0^ 


■ lllitO^ 


,065 


0, 


■ 025it0, 


■ 018 


N 


22 


8 


11 


.927250 


-73, 


.243438 


5, 


.895±0, 


.370 


10, 


.887±1, 


.138 


0, 


■ OOOibO, 


■619 


0^ 


■ 147it0^ 


,038 


0, 


■ 005it0, 


■ 020 


N 


22 


9 


11 


.930068 


-73, 


.312289 


4, 


,740±0, 


,560 


5, 


.598±0, 


.748 


0, 


.255±0, 


■ 267 


0, 


■081it0^ 


,047 


0, 


■ 035it0, 


■ 021 


N 


22 


10 


11 


.930627 


-73, 


.252421 


3, 


.770±0, 


.355 


8, 


.918±0, 


.942 


0, 


■ OOOibO, 


■ 564 


0^ 


■095it0^ 


,051 


0, 


■ 023it0, 


■ 016 


N 


22 


11 


11 


.940693 


-73, 


.289347 


4, 


.091±0, 


.545 


12, 


.358±1, 


.020 


0, 


.650±0, 


■ 237 


0^ 


■956it0^ 


,078 


0, 


■ 463it0, 


■ 028 


N 


22 


12 


11 


.940781 


-73, 


.279369 


4, 


,818±0, 


,457 


13, 


.174±0, 


.957 


0, 


.978±0, 


■311 


0, 


■410it0^ 


,083 


0, 


■ 081it0, 


■ 019 


N 


22 


13 


11 


.940868 


-73, 


.269391 


3, 


,635±0, 


,536 


9, 


.436±1, 


.282 


0, 


.652±0, 


■395 


0, 


■ 212it0^ 


,067 


0, 


■ 056it0, 


■ 024 


N 


22 


14 


11 


.940955 


-73, 


.259413 


4, 


.311±0, 


.433 


9, 


.660±0, 


.727 


1, 


.009±0, 


■ 221 


0^ 


■ 184it0^ 


,054 


0, 


■ 033it0, 


■ 028 


N 


22 


15 


11 


.944087 


-73, 


.298330 


3, 


,994±0, 


,543 


6, 


.894±0, 


.881 


0, 


■ OOOitO, 


■ 212 


0, 


■571it0^ 


,062 


0, 


■ 120it0, 


■ 021 


N 


22 


16 


11 


.944598 


-73, 


.238462 


3, 


,279±0, 


,394 


5, 


.651±0, 


.863 


0, 


■ OOOitO, 


■ 504 


0, 


■ 122it0^ 


,054 


0, 


■ 052it0, 


■ 027 



Note. — All PAH and line measurements are integrated intensities in units of 10~^ W sr~^. The full tabic is published 

in the electronic edition of the journal and includes all of the measurements and uncertainties for the full set of PAH bands and 
the neon lines as well as the values after the "dark settle" correction. 



TABLE 3 
Mean Band Ratios 



Ratio 




N 22 






N 66 






N 76 






N 83 




SW Bar 1 


SW Bar 3 




SMC Bl 




SINGS 




SB 




6.2/11.3 


1, 


■ 165it0, 


■032 


0^ 


■807it0^ 


■037 


1^ 


182it0, 


■ 077 


1, 


■ 153it0, 


■064 


0, 


■882it0.015 


1.240it0, 


■ 028 


0, 


■ 891it0, 


■ 020 


1, 


■ 197it0.276 


1, 


169it0, 


■ 945 


7.7/6.2 


2, 


■ 122it0, 


■071 


2. 


■ 599it0^ 


■ 163 


2. 


,454it0. 


■ 464 


2, 


■ 525it0, 


■ 146 


2, 


■ 525it0.061 


2.428it0, 


■ 089 


2, 


■ 084it0, 


■ 075 


3, 


■ 590it0.638 


4, 


160itl, 


■ 542 


7.7/11.3 


2, 


■487it0, 


■072 


2. 


■ 253it0^ 


■ 121 


2. 


,730it0. 


■ 151 


2, 


■ 968it0, 


■ 118 


2, 


■ 256it0.047 


3.155it0, 


■ 095 


1, 


■ 897it0, 


■ 045 


4, 


■ 176it0.927 


4, 


,684it3. 


■ 979 


8.6/6.2 


0, 


■ 272it0, 


■022 


0. 


■381it0^ 


■022 


0. 


,315it0. 


■ 109 


0, 


■427it0, 


■032 


0, 


■426it0.013 


0.379it0, 


■ 020 


0, 


■ 383it0, 


■ 013 


0, 


■630it0.143 


0, 


,782it0. 


■ 209 


8.6/11.3 


0, 


■ 272it0, 


■034 


0. 


■329it0^ 


■028 


0. 


,369it0. 


■ 064 


0, 


■ 517it0, 


■034 


0, 


■385it0.011 


0.465it0, 


■ 022 


0, 


■ 341it0, 


■ 017 


0, 


■ 723it0.095 


0, 


SOlitO, 


■ 410 


17.0/7.7 


0, 


■ 099it0, 


■004 


0. 


■ 108it0^ 


■009 


0. 


,071it0. 


■ 009 


0, 


■ 063it0, 


■006 


0, 


■ 127it0.003 


0.075it0, 


■ 004 


0, 


■ 123it0, 


■ 014 


0, 


■ 119it0.039 


1, 


,092it4. 


■ 885 


17.0/11.3 


0, 


■ 247it0, 


■007 


0. 


■ 277it0^ 


■017 


0. 


189it0, 


■ 012 


0, 


■ 215it0, 


■017 


0, 


■ 281it0.008 


0.238it0, 


■ 010 


0, 


■ 247it0, 


■ 025 


0, 


■475it0.115 


1, 


,974it6. 


■ 285 


Nelll/Nell 


0, 


■ 541it0, 


■ 107 


5. 


■679it0^ 


■699 


1. 


,905it0. 


■ 260 


0, 


■ 840it0, 


■ 132 


0, 


■394it0.028 


1.043it0, 


■ 233 








5, 


■ 148it4.222 


3, 


148it4, 


■ 800 


6.2/SpAH 


0, 


■ 183it0, 


■004 


0. 


■ 133it0^ 


■008 


0. 


186it0, 


■ 010 


0, 


■ 151it0^ 


■005 


0^ 


■ 145it0.002 


0.169it0, 


■ 004 


0, 


■ 161it0, 


■ 003 


0, 


■ 133it0.019 


0, 


109it0, 


■ 029 




0, 


■ 394it0, 


■006 


0. 


■407it0^ 


■008 


0. 


,398it0. 


■ 024 


0, 


■405it0, 


■008 


0, 


■ 376it0.005 


0.428it0, 


■ 008 


0, 


■ 347it0, 


■ 006 


0, 


■466it0.038 


0, 


,419it0. 


■ 097 


8.6/T,PAH 


0, 


■ 045it0, 


■005 


0. 


■052it0^ 


■003 


0. 


,053it0. 


■ 013 
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correction has a minimal effect on the results. 

To summarize, we find that the 17.0/11.3 ratio is low in 
all SMC regions while the 17.0/7.7 ratio is approximately 
normal compared to the SINGS sample mean. 

3.1.3. Ratios of Individual Bands to the Total PAH 
Emission 

Due to the complexity of comparing among individ- 
ual bands which may vary independently, we attempt 
to simplify the comparison by showing the ratio of each 
band to the total PAH emission. In Figure [T2| we show 
the fraction of the total PAH emission (defined as the 
sum of the 6.2, 7.7, 8.3, 8.6, 11.3, 12.0, 12.6, 13.6 and 
17.0 fim features) carried by each of the major PAH fea- 
tures discussed above. We show the weighted mean of 
the ratio for all of the spectra in a given region where 
the relevant feature is detected at > 5a in both the orig- 
inal and artifact-corrected version of the spectrum. The 
large points show the weighted mean for the original data 
and are connected to a smaller point showing the same 
weighted mean using the artifact-corrected spectra. We 
overplot the mean from the SINGS sample (with an error 
bar showing the standard deviation of the SINGS galax- 
ies) connected with a gray line to guide the eye. We 
can see from this plot that there are deviations from the 
SINGS average: 1) the 7.7 micron complex is weaker, 2) 
the 6.2 and 11.3 fim features are stronger and 3) the 8.6 
and 17.0 ^m bands are weaker compared to the SINGS 
average. We note that the conclusion about the weak- 
ness of the 17.0 /im complex is enforced by the fact that 
it is not a major contributor to the total PAH emission, 
and therefore should be less sensitive to changes in the 
relative proportion of the flux carried by the 6.2, 7.7, 8.6 
and 11.3 /im bands. This agrees with the assertion that 
the weakness of the 7.7 /im feature drives the relatively 
normal 17.0/7.7 ratios seen in the SMC. 

As described in the Introduction, the C-H bands be- 
tween 11—14 fim have the potential to provide informa- 
tion on the structure and hydrogenation of the PAHs. To 
compare the relative strengths of these bands, we make 
a plot very similar to Figure [12] but showing the con- 
tribution of the 11.3, 12.0, 12.6 and 13.6 /xm bands to 
the total C-H band emission (the sum of the emission 
in these features). For the SMC regions we can see from 
this comparison that the 11.3 fim featme carries a similar 
amount of the total C-H band emission as in the SINGS 
galaxies, while the 12.6 iim feature carries slightly less 
and the 12.0 and 13.6 /im features are slightly stronger. 

3.1.4. Summary of the Band Ratio Results 

We have observed clear deviations of the SMC band 
ratios from the average values found in the SINGS sam- 
ple: 

1. For the three brightest PAH bands at 6.2, 7.7 and 
11.3 /im, we find that the 7.7 /im feature is weak 
compared those at 6.2 and 11.3 /im. The 6.2 /im 
band is slightly weak compared to the 11.3 /im 
band. 

2. The 8.6 /im feature appears to be weak relative 
to all of the major bands (6.2, 7.7 and 11.3 /im) 
and a weaker contributor overall to the total PAH 
emission. 



3. The 17.0 /im feature is weak relative to the 6.2 and 
11.3 features. Compared to the SINGS average, 
however, the 17.0/7.7 ratio is only slightly low. The 
17.0 /im feature makes a smaller contribution to the 
total PAH emission in the SMC than in the SINGS 
galaxies. 

3.2. Band Ratios and Radiation Field Hardness 

The ratio of the [Ne HI] line at 15.6 /im to the [Ne 
II] line at 12.8 /im tr aces the hardness of the radiation 
field in H II regions ( [Giveon et al. 2002) and is often 
used in the context of studying PAH emission b ecause of 



the prominence of these two lines in the mid -IR (Madden 
2000| [Smith et aL|2007bl [Gordon et al.|2008[ ). It has been 
suggested that the hardness of the radiation field is the 
driver of the deficit of PAHs observed at low-metallicity 



( |Madden|[2000[ [Gordon et"al]|2008[ ) via enhanced pho- 
todissociation of the molecules. In order to understand 
how the radiation field hardness affects the PAH physi- 
cal state in the SMC, we show the 7.7/11.3 and 17.0/11.3 
ratios as a function of the [Ne III]/[Ne II] ratio in Fig- 
ures [l4| and [TH As will be discussed further in Section [5[ 
these two ratios predominantly trace changes in the ion- 
ization and size distribution of PAHs, respectively. We 
choose these as representative ratios for changes in the 
PAH physical state. We note that none of the other ratios 
discussed here show trends with radiation field hardness. 

From the ratios shown in these figures, we can see that 
there is no clear trend in either the 7.7/11.3 or 17.0/11.3 
with radiation field hardness. The 7.7/11.3 ratio stays 
essentially constant to within ^ 50% over nearly two or- 
ders of magnitude in [Ne III]/[Ne II] ratios. In fact, we 
have found no clear trends of any band ratio with the [Ne 
III]/[Ne II] ratio. The sam e is true for the SINGS galaxies 



with no AGN , as seen by [Smith et al.[ ( [2007b[ ). [Brandl[ 
et al.[ ( [2006[ ) similarly found no trend of the 7.7/11.3 
ratio over an order-of-magnitude in the neon line ratio 
for a sample of starburst galaxies observed with IRS on 
Spitzer. In a resolv ed study of H II regions in M 101, 



Gordon et al. ( 2008 ) found no correlation between their 
ionization index (a combination of the neon ratio and 
the ratio of mid-IR silicon lines) and the 7.7/11.3 ra- 
tio, though they do find that the PAH feature equivalent 
width depends on radiation field hardness. These con- 
sistent results across a wide range of physical conditions 
and metallicity indicate that the radiation field hardness 
does not have a strong effect on the PAH band ratios. 

4. THE PHYSICAL STATE OF SMC PAHS: 
INTERPRETING THE BAND RATIOS 

In Section [3[ we have enumerated the variations and 
average properties of the PAH band ratios we observe 
in the SMC relative to galaxies in the SINGS sample. 
In the following discussion, we use information gathered 
from laboratory and theoretical studies to translate the 
ratios into information about the physical state of the 
PAHs. 

4.1. Interpreting the SMC Results 

Our major conclusions about the band ratios in the 
SMC relative to SINGS are the following: 1) weak 7.7 
/im feature relative to the 6.2 and 11.3 /im features, 2) 
weak 8.6 /im feature relative to all major bands and the 
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total PAH emission and 3) weak 17.0 /im feature relative 
to the total PAH emission. 

We first address conclusion 1), regarding the ratios of 
the 6.2, 7.7 and 11.3 fim bands. As discussed above, the 
6—9 /xm to 11.3 /im band ratios should simultaneously 
trace changes in the PAH ionization and size distribu- 
tion. To aid in interpre ting these rat i os we use a diagram 
based on the model of |Draine fc Li (20011. This model 
utilizes cross sections for ionized and neutral PAHs as a 
function of size based on laboratory and theoretical re- 
sults. The PAHs are exposed to a radiation field with a 
specified strength and spectral shape and their emission 
spectrum is calculated. Figure [TB] shows the band ratios 



from Draine & Li (2001 ) for ionized and neutral PAHs at 
a range of size s exposed to radiatio n fields at 1, 10^ and 



al. (1983) solar neighborhood 



10'^ times the [Mathis et 

radiation field.^ On this plot we show the band ratios 
for the SINGS sample with no AGN and the weighted 
means for the SMC regions. We note that the modeled 
band ratios shown on this plot are for individual PAHs 
of a given size, while PAHs in the ISM should have a 
wide distribution of sizes. Therefor e, in comparing the 



SMC ratios to the models, we use the Draine & Li 



(20011 



results as a guide but not as a way to judge specific PAH 
sizes. 

The relative locations of the SINGS galaxies and the 
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Fig. 16. — Model PAH band ratios from |Drame fc Li| | |2001[ l with the SINGS galaxies without AGN overplotted in black and the weighted 
means of the SMC regions shown in color, 'i'he two sets ot gray Unes show the band ratios for neutral and ionize d PAHs and ea ch at three 
different radiation field strengths (parametrized in units of Ummp i the solar neighborhood radiation field from [Mathis et al.p983[ l . The 
SMC regions fall in the area of this diagram compared to the SINGS galaxies that suggests smaller and more neutral PAHs. 



SMC regions on this plot suggest that SMC PAHs may 
tend to be smaller and more neutral than those in higher 
metallicity galaxies. Within the SMC, N 66 falls closest 
to the SINGS mean, suggesting it has relatively larger 
and more ionized PAHs. SMC Bl, on the other hand, 
has band ratios which put it at the extreme end of the 
diagram, suggesting it has the smallest and most neutral 
PAHs. The spread in the SMC points may be related to 
the progression of the star-formation happening in the 



region: SMC Bl is a quiescen t molecular region (Rubio 
et al.|[l993 Reach et al.||2006 | with abundant molecular 



gas, but no detectable H 11 regions while N 66 has a giant 
H II region with multiple generations of star formati on in 
the vicinity (th ough it also hosts molecular clouds; |Ru- 



bio et al!]|2000 1 . The remaining regions are intermediate 
in the star-tormation properties but generally seem to 
follow a trend where more evolved regions like SW Bar 
3 and N 83 are consistent with slightly larger and more 
ionized PAHs while regions such as SW Bar 1 and N 22 
tend towards smaller and more neutral PAHs. 



It is worthwhile to note that the [Praine fc Li (20011 
model attributes all variation of PAH band strengths to 
either size or ionization. It is very likely that hydrogena- 
tion (in excess of the full aromatic complement of pe- 
ripheral hydrogens) , structure and chemical composition 
also play a role in determining the band ratios. If there 
were changes in structure, we would expect to see vari- 
ations in the relative strengths of the C-H modes. More 
irregular PAHs should show enhanced emission in the 
"duo", "trio" and "quartet" modes as discussed in the 
Introduction. Under some conditions, PAHs can become 
super-hydrogenated — wi th two hydro gen atoms attached 
to a given carbon atom (Le Page et al. 2003). This edge 



structure with two hydrogens per carbon is similar to an 
aliphatic group (— CH2) and is ex pected to produce sim - 
ilarly aliphatic emission features ( Bernstein et al.|[l996 ). 

Our observations in the SMC do not show strong ev- 
idence for changes in PAH structure (towards more ir- 
regular edges) or excess hydrogenation. In Figure 13 we 



show the fraction of the 11—14 /im PAH emission car- 
ried by the 11.3, 12.0, 12.6 and 13.6 /im features. We 
note that in the low-resolution spectra, the 12.6 feature 
is blended with the [Ne II] line at 12.8 /im. For most re- 
gions the 11.3 and 12.0 /xm bands carry a similar fraction 
of the C-H flux as seen in the SINGS sample, however 
the 13.6 /im feature is enhanced and the 12.6 /im feature 
is weaker. Changes in structure that produce more emis- 
sion in "quartet" modes at 13.6 //m would be expected 
to also enhance "trio" modes at 12.6 /im, as both orig- 
inate in irregular edges and extensions. Thus, we can 
not draw a clear conclusion from these band ratios. At 
the signal-to-noise of our spectra, we do not see clear ev- 
idence for any aliphatic features that might result from 
excess hydrogenation. Clear detection of these features 
if they exist and higher signal-to-noise measurements of 
the 13.6 and 14.2 /im bands in the future with JWST will 
enhance our ability to judge whether SMC PAHs show 
changes in structure or hydrogen content. 

We now move on to conclusions 2) and 3) regarding the 
weakness of the 8.6 and 17.0 /xm features in the SMC. As 
discussed in the Introduction both 8.6 and 17.0 fim fea- 
tures are expected to arise in large PAHs. Based on the- 
oretical calculations of PAH spectra, the 8. 6 /im feature 
is thought to trace mainly charged PAHs ( Bauschlicher 
|et al. 2008). There is no clear identification of the charge 
state of the 17.0 /im carrier. The weakness of these two 
features in our spectra provide independent evidence for 
our interpretation that PAHs in the SMC tend to be 
smaller and more neutral than their counterparts in more 
metal rich systems. The particular weakness of the 8.6 
/im feature may be due to the combination of size and 
charge state both favoring less emission in that band. 

4.2. Trends of PAH Physical State with Radiation 
Field Hardness 

In Figures [Ti] and [TS] we show the variation of the band 
ratios as a function of the [Ne III]/[Ne II] ratio. These 
lines arise in ionized gas and their ratio traces the hard- 
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ness of the ionizing radiation field. We find no clear evi- 
dence for variation of the band ratios with radiation field 
hardness over nearly two orders of magnitude in the [Ne 
III]/[Ne II] ratio. PAHs are observed to be depleted in 
ionized gas ( [Povich et "aL]|2007| iLebouteiller et al.|[20ll] ) 
and theoretical studies suggest that the lifetimes of PAHs 
in ionize d gas can be short com pared to the H II region 
lifetimes (Micelotta et al. 2010a). Because of the destruc- 
tion of PAHs in ionized gas, we may not expect the phys- 
ical state of the PAHs to strongly depend on the neon 
line ratio since the PAH emission we observe arises from 
the PDR around the H II region and is thus physically 
separate from the gas traced by the neon lines. 

4.3. Trends of PAH Physical State with Metallicity 

In Figure [17] we show a series of plots illustrating the 
variation of the several PAH band ratios with metallicity. 
In general, the SMC occupies a similar region of the plot 
as other low-metallicity galaxies from the SINGS and 
starburst samples, though there is considerable scatter. 
There is a slight trend towards lower 7.7/11.3 ratios as a 
function of me tallicity, though with many outliers. Smith 
et al. ( 2007b[ ) found a slight trend for galaxies to show 
low 17!0/11.3 ratios at lower metallicity. We show their 
galaxies here, instead using the 17 .Q/Yjpah as a tracer 
of the 17.0 /im feature strength that does not depend 
as heavily on the charge state of the PAHs. The SMC 
17.0 //m fraction of the total PAH emission falls below 
the average for SINGS but above the lowest metallicity 
SINGS galaxy with detected PAH emission (NGC 2915). 

5. DISCUSSION 

5.1. Insights Into the PAH Life-Cycle as a Function of 

Metallicity 

By comparing our measured band ratios with labora- 
tory and theoretical results, we infer that PAHs tend to 
be smaller on average in the SMC than in higher metallic- 
ity galaxies like those represented by the SINGS sample. 
Note that this is not to say that the SMC is devoid of 
large PAHs. We do see a clear 17.0 //m feature in our 
spectra, which requires large PAHs to be present in the 
ISM at some level. The band ratios suggest that the 
proportion of large PAHs is decreased relative to small 
PAHs. If the SMC represents a typical low-metallicity 
galaxy, we argue that the smaller average PAH sizes have 
two major implications for the PAH life-cycle at low- 
metallicity: 1) the change in the size distribution can 
not be the product of processing in the ISM but must in- 
stead reflect changes in PAH formation and 2) the deficit 
of PAHs in low-metallicity galaxies may be related to the 
fact that PAHs are formed with smaller average sizes and 
are therefore more easily destroyed under typical ISM 
conditions. 

5.1.1. Smaller Average PAH Size Due to Formation 
not Processing 

PA Hs can be destroyed i n many ways: photodissocia- 
tion (Allain et al. 1996a|b[ ), interactio ns with energetic 
particles in hot and /or shocked gas ( Micelotta et al. 
2010b|a[) an d interaction with cosmic rays ( Micelotta 
et al. 2011). For all of these processes, our current 
knowledge of PAH destruction suggests that small PAHs 



are destroyed more readily than large PAHs. For in- 
stance, theoretical models of PAH photodissociation in 
the diffuse ISM suggest that PAHs with < 50 C atoms 
are q uickly destroyed under typical interstellar condi- 
tions (jAUain et ;ir]|1996a|b[ |Le Page et al.||2003| while 
large PAHs persist longer iander the same ISM condi- 
tions. Similar trends hol d true for processing by shocks , 
cosmic rays and hot gas CMi celotta et al.|2010b|a] |2011[ ). 
Therefore, since small PAHs are destroyed more easily 
than larger PAHs in the ISM, explaining a deficit of large 
PAHs cannot be the effect of ISM processing unless the 
destruction process turns large PAHs into small PAHs as 
a side-effect 

suggest that this is not the case 



Recent studies by [Micelotta et al. 

The dominant 



2010b I 



lest rue 

tion mechanism for large PAHs is through interaction 
with supernova shocks and they find that "daughter" 
PAHs formed by the fragmentation of the larger grains 
are very quickly destroyed in the shocked gas. Thus, 
since small PAHs are destroyed more easily in all ISM 
processes and they are likely not replenished by the frag- 
mentation of larger PAHs, explaining a deficit of large 
PAHs cannot be the effect of ISM processing. 

We do, however, have ample evidence that processing 
in the ISM can affect the PAH size distribution in the 
sense of selectively removing small PAHs. In the vicin- 
ity of AGN, multiple studies have shown a decrease in 
the 6—9 /^m features relative to the 11.3 um and longer 
wavel ength features ( Smith e t al. "2007b O'Dowd et al. 
20091 [D iamond-St a^^°^^ eke 2010r |Wu et al.| | 2010p . 
Smith et al.| ( ^2007b| ) also found that the 11. 3/17. U ra- 



tio decreases m the presence of an AGN, suggesting a 
stronger contribution from the 17.0 /im feature — a tracer 
of larger PAHs. These trends suggest that the ISM condi- 
tions in the vicinity of the AGN have selectively removed 
the more fragile, smaller PAHs. 
An inte resting counterpoin t to our study of the SMC 



is that by Hunt et al. (2010), who presented a study of 
the PAH band ratios in a sample of blue compact dwarf 
galaxies (BCDs). These BCDs have particularly intense 
radiation fields as well as low metallicities and ver y low 
dust-to-gas ratios (hence decreased dust shielding) 



Hunt 



axies 



et al. (2010) found that the band ratios in these ga 
suggested a size distribution shifted towards larger PAHs 
relative to the galaxies in the SINGS sample (with no 
AGN), particularly in the strength of the 8.6 /xm feature 
and the low ratios of short-to-long wavelength bands. Be- 
cause of the low metallicity of these galaxies, we might 
expect them to show a size distribution shifted towards 
smaller PAHs, as we have seen in the SMC. It may be 
the case in these extreme objects that the radiation field 
is intense enough to remove small PAHs to such a de- 
gree that once again larger PAHs dominate the size dis- 
tribution, regardless of the scarcity of large PAHs. If 
this were the case, we would expect the abundance of 
PAHs relative to dust (i.e. the PAH fraction ) to be very 
low in the BCDs. Indeed, |Hunt et al.[ ([2010[) found that 



the ratio of PAH emission to the total infrared emission 
(EPAH/TIR), which should trace the PAH fraction, is 
depressed by more than an order of magnitude on aver- 
age in these galaxies compared to the SINGS sample. 

If processing in the ISM cannot produce the size distri- 
bution shifted towards smaller PAHs that we observe, the 
only other option is that PAHs are formed with smaller 
sizes in the SMC. Forming PAHs with smaller average 
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Fig. 17. — The band ratios from the 
For the SINGS galaxies we use the metallicities 



SINGS ( [Smith et al.|2007b| l, starburst dErigelbracht et al |200"8ll and SMC samples versus metallicity. 
etallicities tor the central regions from [Moustaitas et al. 1(121)1011 using the calibration of |PiIyugin 



[Th uan {2005J. Th e starburst metallic ities are Tc based and should be consistent with the Filyugin & Thuan (2005) results. They are 
described further in |Engelbracht et al.| ( |2008| ) and do not have uncertainties. We use a metallicity of 12+log(0/H)~ 8.0 for the SMC l |Kurt| 
|fc Dufour|1998^ . Small shitts on the x-axis are introduced for the SMC points so they do not overlap. 



sizes in low-metallicity galaxies could be the result of a 
change with metalHcity in the dominant formation mech- 
anisms (i.e. AGB stars, coagulation or chemistry in dense 
gas — among many other options that have been pro- 
posed) and/or a change in the effectiveness of forming 
large PAHs with the same mechanism. Both of these op- 
tions are plausible, especially given the simple fact that 
the raw material that PAHs form out of is less abundant 
in low-metallicity galaxies, so constructing large PAHs 
may be inherently more difficult. 

5.1.2. Low- Metallicity PAH Deficit as a Consequence 
of Smaller PAH Sizes 

An important consequence of PAHs forming with a dis- 
tribution shifted towards smaller sizes at low-metallicity 
is that a larger fraction of PAHs can be destroyed un- 
der typical ISM conditions. In general, models of PAH 
destruction in the ISM find that b elow a certain PAH 
20 — 50 C atoms, generally; Le Page et al.||200T 



interaction of these two effects could contribute to the 
steepness of the observed drop-off in the PAH fraction 
below a metalhcity of 12-|-log (0/H)~ 8.1 (Engelbracht 
et al.||2005| [Draine et al"]|2007i ). 



5.1.3. Relationship of PAH Emission and Radiation 
Field Hardness 

It has been suggested that the deficiency of PAHs in 
low metallicity galaxies is a product of their destruction 
by harde r and/or more intense radiation fie lds (see for 
instance, |Madden||2000[ IGordon et al.||2008|). We have 



shown that in the SMC the band ratios do not depend 
strongly on the hardness of the ionizing radiation field, 
as one might expect if it were instrumental in destroy- 
ing PAHs. In fact, a wide range of studies have shown 
no correlation between the radiation fiel d hardness and 
the band ratios (Brandl ct al. 2006; Smith et al.||2007b 



size 



Gordon et alpOOS, ,U'Dowd et al.,2009, ,Wu et ar|2010[ ) 



[2007b 
ihyslcal 



2003| PAHs are very quickly destroyed. 'I'he timescale 
for destruction of small PAHs is shorter for essentially 
all destruction mechanisms, as discussed above. With a 
size distribution shifted to smaller sizes, the quickly de- 
stroyed PAHs will represent a larger fraction of the total 
amount of PAHs. Such a situation may naturally lead 
to a PAH deficit in low-metallicity galaxies without the 
need to rely on enhancements in rad iation field hardness 



in line with these two tracers probing distinct pnysical 
regions (H II region vs. PDR). Some studies do show, 
however, a correlation of P AH equivalent width with ra- 



diation field hardness 



2006 Gordon et al. 



equivalent widths in 



Madden] |200q 
2008|, such that PAE 
larder radiation fields 



Beirao et 
have 



al. 

lower 



(Madden 2000 



2008) or delays between 



Gordon et al. 

PAH tormat ion In AGB stars vers us dust formation in 
supernovae (Galliano et al. 2008b I. We note, however. 



various mechanisms can be operating simultaneously- 
i.e. PAHs could be both easier to destroy because of their 
smaller sizes while also being exposed to harder, more 
intense radiation fields in low-metallicity galaxies. The 



This may be 

the product of the correlation between harder and more 
intense radiation fields. Young, massive stars capable of 
creating harder ionizing radiation will also have more in- 
tense radiation fields. The intensity of the radiation field 
may then affect the amount of emission arising in PDRs 
around the HII region. Meanwhile, the characteristics of 
the PAH spectrum are determined by conditions within 
the PDR itself and may be insensitive to the har dness of 
the ionizing radiation (for a similar argument see Brandl 
eral][2006| . 
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5.1.4. Why are SMC PAHs more Neutral? 

In addition to being smaller on average than PAHs in 
the SINGS galaxies, the SMC PAHs also appear to be less 
ionized. One possible reason for lower levels of ionization 
may be that PAHs in the SMC are l argely confined to 
dense regions ( Sandstrom et al^|2010 ). The ratio of the 
ionization to recombination rates should determine the 
PAH charge and is proportional to (Go/rie) x \jTgas- If 
PAHs in the SMC exist in denser regions compared to 
PAHs at higher metallicity, a decreased ionization level 
may be a natural consequence. Another reason why 
PAHs in the SMC may be more neutral is related to 
their smaller average sizes. Small ionized PAHs may 
have a shorter lifetime than small neutral PAHs, since 
ionization in general causes PAHs to be less stable to 
photodissociation. It is possible that PAHs in the SMC 
tend to be more neutral because small, ionized PAHs are 
destroyed even more quickly in the ISM than their neu- 
tral counterparts. More detailed PDR modeling includ- 
ing the PAH ionization state may provide insight into 
why the PAH charge state is different in the SMC than 
in higher metallicity galaxies. 

5.1.5. A Scenario for the PAH Life-Cycle in the SMC 



Sandstrom et al. ( 2010 ) argued that in the SMC the low 
mass traction of dust in PAHs in the diffuse regions com- 
pared to the higher fraction in dense regions suggested 
that PAHs must be forming in dense gas. Their reason- 
ing was that AGB-produced PAHs would primarily be 
deposited in the diffuse ISM which would then condense 
to form dense gas, so enhancements of the PAH fraction 
in dense gas must occur in situ. We propose the following 
(speculative) scenario to jointly explain the abundance 
and physica l state of SMC PAHs inf erred from our study 
and that of Sandstrom et al. (20101: 1) PAHs are form- 
ing in regions of dense gas in the SMC; 2) these PAHs 
are formed with a smaller average size possibly due to 
the lower abundance of carbon available in the ISM; 3) 
as they emerge from dense regions into the diffuse ISM, 
the PAHs interact with UV photons, shocks and cosmic 
rays, and a fraction of the PAHs are preferentially de- 
stroyed on the small end of the PAH size distribution. In 
this scenario, a large fraction of the PAHs in the SMC 
are destroyed in the diffuse ISM, causing a PAH deficit, 
not necessarily because of harsher conditions, but be- 
cause the average PAH starts off smaller and thus more 
"fragile" in the SMC than the average PAH in a higher 
metallicity galaxy. 

6. SUMMARY & CONCLUSIONS 

We have presented a study of the PAH physical state 
diagnosed through the ratios of the 6—17 fim bands in 
several regions of the SMC. In comparing our measured 
band ratios to the SINGS sample, excluding galaxies with 
AGN, we find several distinct trends: SMC PAHs show 
a weak 7.7 /im feature relative to the other dominant 
PAH bands at 6.2 and 11.3 /im and the features at 8.6 
and 17.0 fim are weaker contributors to the overall PAH 
emission. We find no trend of the band ratios as a func- 



tion of radiation field hardness as traced by the mid-IR 
[Ne III]/[Ne II] ratio. 

Interpreting the band ratios in light of laboratory and 
theoretical results on the PAH spectrum, we suggest that 
SMC PAHs are smaller and more neutral than their coun- 
terparts in higher metallicity galaxies. Several lines of 
evidence argue for this interpretation: the 7.7/11.3 and 
7.7/6.2 ratio suggests that PAHs ar e smaller and more 
neutral based on comparison with the Draine & Li ( 2001 1 
model and various theoretical and laboratory studies; the 
8.6 /im band, which is thought to be dominantly emitted 
by large, charged PAHs is particularly weak in the SMC; 
and the 17.0 /xm feature — which is thought to arise in 
skeletal vibrations of large PAHs — is also weak relative 
to the total PAH emission in the SMC. 

We argue that a size distribution shifted towards 
smaller average sizes can not be the product of ISM 
processing by UV photons, shocks, cosmic rays or hot 
gas, since all of those processes more efficiently destroy 
small PAHs. A variety of literature results suggest that 
ISM processing creates a size distribution shifted towards 
larger PAH sizes, the opposite trend we see here. If ISM 
processing can not create the size distribution we ob- 
serve, the likely explanation is a change in how PAHs 
form. This could be a change in the mechanism of PAH 
formation at low metallicity or a change in the efficiency 
of creating large PAHs. 

Our observations imply that the PAH deficiency at 
low-metallicities may not be due to enhancements in the 
hardness or intensity of radiation fields, but rather that 
PAHs at low-metallicity start off inherently more "frag- 
ile" because of their small sizes so a larger proportion of 
them can be destroyed under typical ISM conditions. 
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Fig. Al. — Several spectra from the map of SW Bar 1 illustrating various degrees of mis-match at the SL2/SL1 and LL2/LL1 overlap 
regions. The uncorrected orders are shown in color — SL2 in orange, the SLl in purple, the LL2 in green and the LLl order in gray. The 
offsets in the overlap regions are due to a time- and positionally varying dark level on the detector. The black spectrum overlayed shows 
these same regions after the correction described in this Appendix. The correction derived from the inter-order regions greatly improves 
the mis-match at the SL2/SL1 overlap, but it may oversubtract real 7.7 fim PAH emission. 



APPENDIX 

ARTIFACTS IN THE IRS SPECTRAL MAPPING AND THEIR MITICATION 

Description of the Artifact 

In reducing the S'*MC observations, we found that the spectra were affected by an artifact caused by a time- and 
positionaUy variable background level on the detector. This artifact occurred in both SL and LL observations. It 
appears to be similar to the "dark settle" artifact that can be important for high-resolution observations with IRS. 
The artifact causes mis-match in the regions where the orders overlap and was first identified because it causes the 
long wavelength end of SL2 and the short wavelength end of SLl to be offset from each other, sometimes significantly. 
Since this wavelength region covers the 7.7 /im PAH feature, understanding the effects of this artifact was crucial for 
our study. In Figure |Al| we show several spectra from the SW Bar 1 region illustrating the mismatches at SL and LL 
caused by the artifact. Similar effects are seen in all maps, except for the deep SMC Bl #1 SL map, which used 60 
second ramps and does not appear to suffer from this artifact. 

We found that the shape and magnitude of the offsets can generally be explained by residual background emission 
that is apparent in the inter-order region of the BCDs. The residuals persists after removal of the "off" observations 
and are much larger in magnitude than the difference in zodiacal light or MW foreground between the map position 
and the "off" position. The average properties of this dark level can be seen in the median background-subtracted 
BCD for each mapping AOR — created by subtracting the median background image (created in CUBISM) from each 
of the BCDs in a mapping AOR and then combining the stack of all of the BCDs to generate one median background- 



subtracted BCD. Figure X2 shows these for each of our SL mapping observations and Figure [Xs] shows these for each 
of our LL mapping observations. In these Figures we have masked out the orders and the peak-up arrays to show the 
residual background more clearly. 



In addition to the spatial structure of the background, there is also some degree of time- variability. In Figure X4 



we illustrate the time variations of the residual background by displaying the median of a ~5 pixel wide strip of the 
background subtracted BCDs versus BCD number (essentially a time sequence of the observations). 

Correction of the Artifact Derived from Inter-order Background Level 

To correct the effects of the artifact, we use the inter-order regions to determine the spatial shape and time- variability 
of the residual background. This approach is very similar to that done by Starkey et al. (2011, in prep). Because 
of pixel-based noise, we generally need to average BCDs together as a function of time to improve our determination 
of the background. At SL wavelengths the presence of the "peak up" arrays on the right-side of the BCD make the 
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Fig. A2.— The median of all background subtracted BCDs for the SL mapping AORs in a) N22 North, b) N66 North, c) N 76 #1, d) 
N 83, e) SW Bar 1 and f) SW Bar 3. The SLl, SL2 and SL3 orders and the peak-up arrays have been masked to more clearly show the 
background structure. These images highlight the non-time variable characteristics of the residual background, which are similar in each 
mapping AOR. These include a higher background level at the top of the detector (~8 fj.m in the SLl order and SL3) and in the middle 
of the detector (at ^ 5-6 fim in the SL2 and ~ 12 /im in SLl). 




-4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 
BCD Units 



Fig. A3.— The median of aU background subtracted BCDs for the LL mapping AORs in a) N22, b) N66, c) N76, d) N83, e) SW Bar 1 
and f) SW Bar 3. The LLl, LL2 and LL3 orders have been masked out to show the residual background more clearly. These images 
highlight the non-time variable characteristics of the residual background. For the LL maps, the residual background is concentrate d in 
the bottom right of the detector. This causes an offset at the long wavelength end of LL2, which can be seen for SW Bar 1 in Figure [Al] 
when the LL2 and LLl orders do not agree in their overlap region. 
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Fig. A4. — The median value of the background subtracted BCD in a ~5 pixel wide box (overlayed in red on the median BCD on the 
left) as a function of BCD number, which generally corresponds to time in the AOR. Panel a) shows the N66 South SL mapping AOR that 
includes 600 BCDs and panel b) shows the N 66 LL mapping AOR that includes 686 BCDs. The images have been smoothed with a 3x3 
pixel Gaussian to suppress noise. These images highlight the time-variable characteristics of the residual background level. 

available inter-order regions very small, since scattered light from the peak up arrays affects the nearby inter-order 
regions. We are therefore limited to the small strip between the SLl and SL2/S L3 o rders to characterize the vertical 
spatial dependence of the residual background (this region is highlighted in Figure [A4| . For the SL orders we determine 
the median BCD value as a function of time and vertical row in a box ~5 pixels wide by 11 pixels high by 3 BCDs 
long. We then create a correction image for each BCD where all columns contains this median. Finally, we apply the 
SL flatfield to get the appropriate value within the orders. 

At LL wavelengths, the lack of peak up arrays make this process much simpler. In addition, there is far less time 
variability in the LL residual background so we use a simpler correction derived from the median BCDs shown in 
Figure 1^3] We mask the LL orders and fit a horizontal polynomial to each row of the BCD, filling in each row of the 
correction image with this fit. We then apply a vertical polynomial smoothing in a 10 pixel box and finally use the LL 
flatfield image to obtain the correction in the LL orders. For both the SL and LL we subtract these correction images 
from the original BCDs. 

The effect of the correction on several spectra can be seen in Figure [XT] in black. The agreement between orders in 
the overlap region is essentially always improved after the correction. Most of this correction comes from removing 
flux in the 7—8 /im region, coincident with the 7.7 PAH feature. Because of the limited information about the spatial 
structure of the residual background and the large effect it has on the PAH feature, we use both the corrected and 
uncorrected spectra to analyze the PAH band ratios in this paper. Future work my reveal the source of this residual 
background and improve our determination of the PAH feature strengths. 

To conclude we make a few remarks about the sort of observations this artifact is likely to affect. Generally, the 
residual background level is not an issue in typical IRS observations when either (1) an in-slit background can be 
subtracted, as in the case of point source observations; (2) the background observation is very close in time to the 
mapping observation, for example when a "nod" position can be used for background subtraction; or (3) when the 
object being mapped is bright enough that these low-level variations are negligible. In the SMC we are observing faint 
emission using the full IRS slit and our "off" observation occurs after ^600 pointings. Thus, the S^MC observations 
present a somewhat unique situation where the residual background artifacts are important to remove. 

Facilities: Spitzer () 
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